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A STUDY OF CRACK CLOSURE IN FATIGUE 

t>y 

T. T. Shih and R. P. Wei 
LEHIGH UNIVERSITY 
Bethlehem, Pennsylvania 

SUMMARY 


Crack closure phenomenon in fatigue was studied by using a Ti-6A1-4V titanium 
alloy. The occurrence of crack closure was directly measured by an electrical-potential 
method, and indirectly by load-strain measurement. The experimental results showed 
that the onset of crack closure depends on both the stress ratio, R, and the maximum 
stress intensity factor, K max . No crack closure was observed for stress ratio, R, 
greater than 0. 3 in this alloy. 

A two-dimensional elastic model, was used to explain the behavior of the recorded 
load-strain curves. Closure force was estimated by using this model. Yield level 
stress was found near the crack tip. Based on this estimated closure force, the crack 
opening displacement was calculated. This result showed that onset of crack closure 
detected by electrical-potential measurement and crack-opening-displacement measure- 
ment is the same. 

The implications of crack closure on fatigue crack are considered. The experi- 
mental results show that crack closure cannot fully account for the effect of stress 
ratio, R, on crack growth, and that it cannot be regarded as the sole cause for delay. 
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I. INTRODUCTION 


The importance of fatigue in determining the serviceable life of engi- 
neering structures has been well recognized [1], In modern high-performance 
structure designed for finite service life, fatigue crack growth occurs over a 
significant portion of the useful life of the structure. Information on the 
kinetics of fatigue crack growth, therefore, becomes an important aspect of 
material selection and design. 

It has been shown that the crack tip stress intensity factor, K, defined 
by linear elasticity, is the most appropriate parameter for characterizing 
crack growth [2,3,4], Several empirical relationships for correlating fatigue 
crack growth data have been suggested [2, 5-8]. Paris et al [7] and Paris and 
Erdogan [8] suggested that the primary variable of importance for fatigue crack 
growth is the range of the crack tip stress intensity factor, AK, and that the 
rate of fatigue crack growth, A a/ A N, may be described by a relationship given 
by Equation 1. 

^ = A <AK) n (1) 

A and n are empirical constants. The mean value, or the maximum value of K 
(•^mean or -^max) ™ a £i ven l oa d cycle was thought to have only a secondary 
effect on crack growth. Forman et al [6] recognized that the maximum value 
of K, K max , can have a significant effect on fatigue crack growth, particularly 
at maximum K levels approaching (the critical stress intensity for failure). 
They suggested that Equation 1 be modified so that the rate of fatigue crack 
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growth would approach infinity as approaches K c , Equation 2. 

Aa = A ( AK) n 
AN (l-R)!^ - AK 

The stress ratio, R, is the ratio between the minimum and maximum K levels 

in a load cycle; R = Reasonable correlation with fatigue crack 

growth data in the range of 10“® to 10”® inch per cycle (2. 54 x 10”® to 2. 54 x 

10” 2 mm per cycle) has been demonstrated [9]. 

The implicit assumptions in both Equations 1 and 2 are (a) that only the 
tensile portion of the load cycle is effective in producing fatigue crack growth, 
and (b) that there is no crack closure over the tensile portion of the load cycle, 
so that the full tensile range of loading is effective. In a recent series of 
experiments on a 2024-T3 aluminum alloy, Elber observed that the load versus 
crack-opening-displacement curves exhibited a nonlinear region at the lower 
load levels, as indicated schematically in Figure 1 [10], This behavior was 
interpreted in terms of crack closure, that is, physical contact between the 
fracture surfaces produced by fatigue. Elber suggested that the crack is 
closed at the tip over the lower porition of the loading cycles and becomes 
open only after the applied stress exceeded a level s op (Figure 1) and that 
fatigue crack growth can occur only during that portion of the loading cycle in 
which the crack is fully open. Based on this suggestion, an effective stress 
range, A S e £j, and an effective stress range ratio, U, were defined. 

A s eff _ s max - s op 
U = AS " S max - S min 

S max and S m j n are the maximum and minimum values of the applied stress in 
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a given cycle; and S Q p is the crack opening stress. U can be defined equiva- 
lently in terms of the effective stress intensity range, AK e ff, and Ak. 

AKeff _ K ma x - K 0 p 

K _ A T/ ~ K • ^ 

A max ^min 

K Q p is the crack opening stress intensity factor corresponding to S Q p. Elber 
further suggested that the crack growth relationship be written in the following 
form: 

2^- = C (AK eff ) n = C (UAK) n (5) 

Based on a limited range of data, Elber suggested that the effective 
stress range ratio, U, is a function of the stress ratio, R, and is independent 
of ®max> or Kmax* For the 2024-T3 aluminum alloy sheet, U is given simply 
by Equation 6, for R value ranging from -0. 1 to 0. 7. 

U = 0.5 + 0.4R -0.1$ R <0.7 (6) 

Using this empirical result, Elber showed that Equation 5 provided a better fit 
to the experimental data than either Equation 1 or Equation 2 [10], Because 
Equation 6 was based on a very limited amount of data and because much of the 
data was obtained under conditions of net-section yielding, this relationship for U 
should be considered only as a preliminary estimate and must be verified by 
additional experimentation. Thus far, no independent verification of this result 
has been reported, although some direct evidence for closure has been given by 
Buck et al [11-13]. Nevertheless, Elber's concept of effective stress range has 
been adopted by several authors [14,15], and the expression for U (Equation 6) 
has been accepted and used without question by others [16-18], 
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The principal purpose of this work is: (a) to provide independent verifi- 
cation of the crackclosure phenomenon; (b) to determine the effects of stress 
ratio, R, and maximum stress intensity factor, K max , on closure; and (c) to 
develop an analytical model for estimating the magnitude and distribution of 
the closure forces. The experimental and analytical studies are carried out 
within the framework of linear-elastic fracture mechanics, A Ti-6A1-4V alloy 
sheet is used in the experimental studies on closure. Crack closure is mea- 
sured directly by means of an electrical potential technique, and indirectly by 
measurements of strains in the neighborhood of the crack tip. The implications 
of this study to the understanding of several practical problems in fatigue crack 
growth are considered. 

Acknowledgment is extended to Mr. J. H. FitzGerald for his invaluable 
contributions to the experimental work. 
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II. MATERIAL AND EXPERIMENTAL WORK 


A. Material and Specimen 

A 0. 2- inch- thick (5.08 mm) mill annealed Ti-6A1-4V alloy plate was used 
in this investigation. The chemical composition, and longitudinal and trans- 
verse tensile properties of this alloy are given in Table I. 

Three- inch- wide (76.2 mm) by 16-inch-long (406 mm) center-cracked 
specimens (Figure 2), oriented in the long transverse (TL) direction, were 
used in the fatigue crack growth and crack closure studies. The initial center 
notch, about 0.4 inch (about 10 mm) long, was introduced by electro-discharge 
machining (EDM). The specimens were precracked in fatigue either at leads to be 
used in the subsequent studies, or through a decreasing sequence of loads that 
terminated at these loads. The precracking procedure provided fatigue cracks 
about 0. 08 inch (about 2mm) in length from the ends of the starter notches cor- 
responding to a half-crack length of about 0. 3 inch (7. 6 mm), and ensured that 
the subsequent tests will be unaffected by the starter notch geometry and by the 
residual stress produced by EDM. Both precracking, and the crack growth and, 
crack closure experiments were carried out in a 100,000-lb. capacity MTS 
system closed-loop electrohydraulic testing machine. Load control was esti- 
mated to be better than + 1 percent. 

The stress intensity factor, K, for the center-cracked specimen was 
computed from Equation 7: 

K = Sn a sec (it a/W) (7) 

BW 

where P = applied load, B = specimen thickness, W = gross specimen width, 
and a = half crack length. The secant term is a correction for finite specimen 
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width [19] which closely approximates the series correction obtained by Isida 

[ 20 ]. 


B. Crack Monitoring System 

An electrical potential method was used for monitoring crack growth and 
closure utilizing a continuous recording system. This method is based on an 
increase in electrical resistance of the specimen with crack growth. A con- 
stant current is applied to the specimen, and changes in electrical potential 
are measured between fixed points above and below the crack, Figure 3. A 
schematic diagram of the system is shown in Figure 4. The relationship 
between the electrical potential, V, half crack length, a, and the distance 
between the two fixed measurement points, 2y, for the central cracked speci- 
men is given by Equation 8: 


V(a) 

V(ao) 


cosh -1 ( 

'cosh TTy/Wl 

icos 7 ia/W / 

cosh -1 j 

[cosh TTy/Wi 

^cos 7 ra 0 /W/ 


( 8 ) 


This equation, derived by Johnson [2i] was used as an analytical calibration 
curve. The electrical potential method provides measurements of average 
crack length through the thickness, as opposed to the visual methods which 
give measurements of the crack length at the specimen surface only. This 
method has been shown to agree well with other crack measuring techniques 
for a number of materials tested in various environments, provided that crack 
tunneling is taken into account [22], 

This method of crack measurement has several advantages. It permits 
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measurements of crack length while the crack is completely covered, thus giving 
complete freedom for using environmental chambers which may completely 
cover the crack area. Since changes in electrical potential reflect electrical 
shorting, vis-a-vis, physical contact, across the crack surfaces, these changes 
provide direct measures of crack closure. The other important advantages are 
that it gives a continuous measurement of crack length as a function of time and 
that it permits continuous recording of load-potential curves. By using a 
working current of about 2. 1 amperes, this system provided an average mea- 
surement sensitivity of about 0. 0025 inch (0. 064 mm) in half crack length, a, 
per microvolt (juv) changes in potential (that is, 0. 0025 inch/yuv, or 0. 064 mm 
juv) for these specimens. Resolution is better than 0. 001 inch, or 0. 025 mm. 

For autographic recording of the load versus potential data used in the 
crack closure studies, an alternate amplification procedure was used in place 
of the system shown in Figure 4. The electrical potential signal from the 
specimen was applied directly to the input of a d-c amplifier in the MTS 
Systems testing machine. The output from this amplifier was reduced by a 
preset d-c signal, and the difference signal was amplified further in the x-y 
recorder and recorded. Because these amplifiers were designed for signal 
levels in the millivolt range, and the changes in electrical potential were in 
the microvolt range, the background noise tended to be higher with this 
method, The reduced signal to noise ratio, however, did not seriously affect 
the closure results (see Figure 6, for example) and was acceptable in view of 
the improved convenience in measurement. 
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C. Environment Control System 


Preliminary experiments indicated that an insulating oxide layer formed 
on the fracture surfaces of specimens tested in air. This insulating oxide 
layer interferred with the measurement of crack closure by the electrical 
potential method, and gave values that underestimated the extent of closure. 

To circumvent this problem, all electrical potential measurements of closure 
were carried out on specimens tested in dehumidified argon that had been further 
purified by a titanium sublimation pump. A schematic diagram of the environ- 
ment control system used in these experiments is shown in Figure 5. 

Purification was achieved by passing ultrahigh purity argon through a 
molecular sieve drier, and a series of cold traps at about -140°C. Further 
purification was obtained by passing the argon through a titanium sublimation 
pump (TSP) in line with the gas purification system. The TSP was operated 
as a getter, and was connected directly to the environment chamber through a 
high conductance tube. To reduce possible back diffusion of impurities, the 
argon was discharged through an additional cold trap and a silicone fluid trap. 
The effectiveness of this purification system has been demonstrated by Wei 
and Ritter [23]. 

D. Experimental Work 

Two independent methods were used to study the crack closure phenome- 
non. Crack closure was measured directly by using the electrical potential 
method. Indirect measurements of crack closure were made by measuring 
strains in regions near the crack tip. To simplify experimentation, only one of 
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the methods was used on each test specimen. The effects of both R and 
were examined. 

Closure experiments were carried out, at suitable crack length intervals, 

on specimens that have been fatigued under constant load- amplitude cycling at 

selected stress ratios, R. Stress ratios of 0. 05, 0. 1, 0. 2, 0. 3 and 0. 5 were 

used with the electrical potential method. For the strain measurements, only 

R of 0. 05 was utilized. A constant maximum load, P max » of 12, 000 lb 

(53. 38 kN) was used for all the tests except for those at R = 0. 05. For the 

tests at R = 0. 05, P max 10, 000 lb. (44.48 kN) was used. Fatigue cracks 

were extended from an initial length of about 0. 3 inch to about 0. 9 inch (or, 

from about 7.6 to 23 mm). The corresponding values of ranged from 

about 15 to 33 ksi-in 5 (16. 5 to 36. 3 MN-m - ^/^) f or the test at a P max of 

10, 000 lb. (44.48 kN), and from about 17 to 40 ksi-in 5 (18. 7 to 44 MN-m - ^/^) 

for tests at P„„. of 12, 000 lb. (53. 38 kN). 
max 

1. Electrical Potential Measurements 

Electrical potential measurements of crack closure were carried 
out on specimens that were fatigued in dehumidified argon (see section on 
Environment Control System). Fatigue cracks were extended to various selec- 
ted lengths under constant load-amplitude cycling at 5 Hz. for prescribed R 
and P m ax* Cyclic loading was interrupted at selected crack length intervals, 
with the load at P max > and the closure experiments were performed. The fol- 
lowing sequences were utilized and constitute one set of closure experiments: 
(a) the specimen was unloaded from P max to P m jn, and then reloaded to P m ax 


12 



(b) the specimen was then unloaded from P max into compression to about 4, 000 
lb. (17.79 kN compression), and reloaded to P max ; and (c) finally, step (a) was 
repeated. These procedures were carried out using the one-cycle loading feature 
on the MTS machine operated at 0. 01 Hz. Autographic recordings of load 
versus change in electrical potential (vis-a-vis, change in apparent crack 
length - crack closure) were made during each of the aforementioned steps. 

2. Strain Measurements 

For convenience, specimens for the strain measurements were 
fatigued in air. Metal foil strain gages, with 0. 015 inch (0. 38 mm) gage length, 
were used. Up to 4 strain gages were mounted at various locations on the 
specimen surface adjacent to the line of intended crack prolongation (see 
Table 2 for specific locations for the strain gages). Fatigue cracks were 
extended under constant amplitude cyclic loading at 5 Hz. , with P max = 

10, 000 lb. (44. 48 kN) and R = 0. 05. Cyclic loading was again interrupted at 
selected intervals (to approximate those of the electrical potential measure- 
ments) for the closure studies. Here, only the load sequence P max to 4, 000 lb. 
compression (-17. 79 kN) back to P max was used. Autographic recordings of 
load versus strain were made during these unloading and reloading cycles. 
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III. RESULTS AND DISCUSSIONS 


In this section, the experimental results will be summarized. The 
results of the electrical potential and strain measurements will be considered 
separately. The effects of stress ratio, R, and maximum stress intensity 
factor, K max , on crack closure will be discussed. Detailed interpretations of 
these experimental results and correlations with an approximate analysis will 
be given in the following section on analytical modeling. 

A. Electrical Potential Measurements 

Typical load versus change in electrical potential curves are shown in 
Figure 6-8 for the sets of crack closure experiments described in paragraph 
D-l of Section II. Figure 6 depicts those for loading sequence (a) — P m ax to 
p min t0 p max! Figure 7, those for sequence (b) — Pmax to -4000 lb. (-17. 79 kN) 
to P m ax; and Figure 8, those for sequence (c) — P m ax to P m in to p max — 
following the sequence into compression. Each of the curves corresponds to 
a different crack length, and thus corresponds to a different value of K max . 
Sensitivity ranged from about 0. 005 in/pv. (0. 13 m »/pv. ) to 0. 002 in/pv 
(0. 05 mm/pv), corresponding to crack lengths, a, of 0. 3 to 0. 8 in. , respectively. 

Changes in electrical potential reflects electrical (and, hence physical) 
contact between the fracture surfaces, and thus provide a direct measure of 
crack closure. For a crack in an elastic medium, the crack surfaces are expected 
to be completely separated (open) under an externally applied tensile load, and to 
be incomplete contact (fully closed) in compression. Load versus change in electrical 
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potential curves for this idealized case are expected to follow the behavior 
indicated by Figure 9. In the tensile region, the electrical potential assumes 
some value V(a) corresponding to the prevailing crack length. As the applied 
load is reduced through zero into compression, the electrical potential under- 
goes a step-wise change from V(a), and assumes a value V(a 0 ), corresponding 
to that for some initial finite-width notch of length a 0 , or a value corresponding 
to the uncracked specimen. In reality, if crack closure occurs, it is expected 
to proceed from the crack tip and extend gradually back towards the initial 
notch. This gradual change is reflected in the actual experimental data, 

Figures 6-8. Initial deviation from V(a) can be identified with the onset of 
crack closure. This point can be identified with the crack opening stress, 

S 0 p, used by Elber [10], and will be utilized in computing U. (Note that, aside 
from the identification with the onset of crack closure, no further physical 
significance is assumed or implied. Further discussion of this point will be 
given in a later section. ) The influences of R and K max on the onset of crack 
closure will be considered separately in a later section. 

It is useful to consider the nature of the load versus change in electrical 
potential curves in some detail, which serves to qualify the experimental 
results and provides some physical insight into the crack closure phenomenon. 
It is important to recognize that the electrical potential measurements provide 
measures of the average (through the thickness average) crack lengths, and 
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that the identification of a change in electrical potential with an increment of 
crack closure is based on the assumption that the electrical conductance across 
the fracture faces approximates that of the undamaged material. It is believed 
that this assumption was reasonably well met, and was correct in the case of 
compression loading based on the agreement between the electrical potential 
under compression loading and that of the uncracked specimen. Irrespective 
of this assumption, the point for the onset of crack closure can be identified 
with good accuracy. 

From the experimental data, it can be seen that the paths of the unloading 
and reloading curves were somewhat different. The differences are particularly 
apparent in those cases where the specimens were loaded into compression. 

The differences were caused in part by short-term drift and electrical noise in 
the measurement circuit. The principal difference, Figure 7, was produced, 
most likely by crushing of the fracture surfaces, and by refracturing of regions 
of the fracture surfaces that had become "cold welded" during compression. 

B. Strain Measurements 

For an idealized crack in an elastic medium, as considered in the pre- 
vious section, the load versus strain curves will exhibit two types of behavior 
depending on whether the strain gage is located ahead of, or behind the crack 
tip. For a gage that is located ahead of the crack tip, the idealized behavior 
would be that shown in Figure 10(a). In the tension region, the behavior is 
linear and reflects the stress (strain) concentration effect of the crack, hi 
compression, the slope of the load-strain curve corresponds to that of an 
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uncracked specimen. The abrupt change in slope at zero load reflects the 
abrupt change in stiffness as this idealized crack changes from a fully open to 
a fully closed configuration. Figure 10(b) illustrates the idealized load-strain 
behavior, when the strain gage is located behind the crack tip. In the com- 
pression region, the load strain curve again corresponds to that of an uncracked 
specimen. In the tension region, since the strain (stress) is essentially unloaded 
by the presence of the crack, changes in strain with load become very small. 
Similar to the case represented by Figure 10(a), the abrupt change in slope at 
zero load is again that associated with opening and closing of the crack. 

The load versus strain behavior for cracks in real material is repre- 
sented by the curves shown in Figures 11. The load was again cycled from 
P max into compression, and then back to P max « The curves represent mea- 
surements made on two separate gages, and reflect changes in load-strain 
behavior as the crack tip approached and then moved away from the gage. The 
substantial deviations in behavior from those of the idealized crack again reflects 
the gradual nature of the closure process. The onset of closure can be detected 
most easily when the crack tip is close to the strain gage, and is determined 
by the deviation from linear behavior in the tension region. The loads at which 
the onset of closure began correlate well with the loads obtained from the 
electrical potential measurements. A more detailed discussion of these results 
is best made in terms of an analytical model , and will be deferred. 
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C. The Effects of Stress Ratio and K max on Crack Closure 

On the basis of the electrical potential data, Figures 6 and 7, the effect 
of stress ratio, R, and the maximum stress intensity, K max , on crack closure 
may be considered. Crack closure is considered to begin when the electrical 
potential deviates from the value at full load. The point for the onset of crack 
closure is consistent with that defined by Elber [10] on the basis of crack opening 
displacement measurements. The results are shown in Figure 12, in terms of U 
(the effective stress range ratio, U, as defined by Elber). It is to be emphasized 
again that U is to be interpreted here only in terms of the onset of crack closure. 
No further physical significance is assumed or implied. 

U and P Q p, corresponding to the onset of crack closure, are shown in 
Figures 12a and 12b as functions of K max for various R. These results clearly 
indicate that the onset of closure is a function of stress ratio, as well as a func- 
tion of K max and hence of crack length. For R greater than 0. 3, no crack closure 
was observed (i. e. P Q p < Pjjj^), and U = 1. For values of R below 0.3, U is a 
function of K max above certain K max values, and is independent of K max (i. e. U 
approaches 1) below these K max values, Figure 12a. It is dependent on K ma x in 
an intermediate range, and becomes mildly dependent on K max at high K max 
levels. The latter region corresponds to the range in which Pop becomes nearly 
constant, Figure 12b. These results are qualitatively consistent with the physical 
model (see next section) that suggested an influence of stress ratio, R, and are 
consistent with the probable dominance of the surface plastic zones on closure. 
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The results are, however, in contrast with those reported 
by Elber [10] for a 2024- T3 aluminum alloy. The differences, of course, can 
be attributed in part to the differences in materials. The principal causes, 
however, reside in Elber 's data. First, because of the limited range of data 
(obtained mostly at high K max levels) and of the inherent uncertainties in the 
results, Elber's assertion that U is independent of K max cannot be regarded as 
being conclusive. Second, most of Elber’s data were obtained 
under conditions above net section yielding. As such, the utility of the data and 
of the suggested relationship between U and R is questionable. Further independ- 
ent verifications are needed. 

D. Discussion 

Based on the experimental observations, the following processes for 
crack closure are envisioned. The physical basis for closure can best be des- 
cribed by comparison with idealized cases. For an idealized stationary crack 
(infinitesimally thin slit) in an elastic medium, the crack surfaces are expected 
to be completely separated (open) under an externally applied tensile load. 

The separation between the crack faces (crack opening displacement), can be 
defined by linear elastic analysis [24]. Upon unloading to zero load, the crack 
opening displacement would return to zero along the full length of the crack. 

For an idealized stationary crack in an elastic-plastic medium, the crack 
opening displacements (crack contour) under tension are expected to be larger 
than those of the elastic case at the same load, as a result of plastic 
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deformation ahead of the crack tip [25]. Because of this plastic deformation, 
residual crack opening displacements will remain following the reduction to zero 
load, and no crack closure would be indicated. 

Real cracks, however, are formed as a result of highly localized defor- 
mation and separation at the crack tip, and are expected to behave quite dif- 
ferently. It is envisioned that the process of deformation and separation (for 
example, that suggested by Krafft [26, 27]) followed by compaction, would pro- 
duce a highly deformed layer. This layer, defined as "a layer of residual ten- 
sile strain left in the wake of the crack tip, " [10] delineates an effective crack 
contour with crack opening displacements that are, in all probability, smaller 
than those of an equivalent elastic crack. On unloading, the surfaces of this 
layer may come into contact and give rise to the observed results. Of course, 
the degree of compaction depends on the ratio between the minimum and maxi- 
mum loads in a cycle, that is, on stress ratio, R. Hence the onset of crack 
closure is expected to be a function of R. Since the extent of the residual defor- 
mation is dependent on K max , the onset of crack closure is expected to be a 
function of K max , also. The observations of experimental results support this 
model. 
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IV. MODELING FOR CRACK CLOSURE 


A two-dimensional elastic model is used to assist in the interpretation of 
experimental results. It is recognized that such a model cannot truly repre- 
sent the real physical problem, which is one that involves elastic-plastic 
behavior, with unloading, and is most probably three-dimensional in nature 
[28], Nevertheless, it is felt that some semi-quantitative understanding can be 
obtained. In the following sections a brief description of the model is given, 
and the qualitative features of the model are discussed in terms of the observed 
load-strain behavior. On the basis of the model and the experimental data, 
estimates of the crack closure forces and the load versus crack opening displace- 
ment curves are made and discussed. 

A. Analytical Model 

The model chosen is that of a notch-like crack in a homogeneous, iso- 
tropic elastic body loaded by remote tension and by distributed compressive forces, 
representing forces produced by crack closure, on the crack faces as shown in 
Figure 13. The magnitude and distribution of the distributed forces are expected to 
depend on the degree of closure. It is assumed that the crack faces remain 
separated, and the only manifestation of closure is the artifical introduction of 
closure forces. It is further assumed that the closure forces are uniformly 
distributed through the thicxness (that is, in the z-direction), and vary only in 
the x-direction, or along the crack length direction. The problem is then 
treated as one of generalized plane stress. Based on these assumptions, the 
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stresses (strains) at any point can be obtained by superposition, Figure 13. 

This model is analogous to the strip yield model proposed by Dugdale [29], and 
is expected to provide reasonable approximate solutions for stresses (strains) 
away from the crack tip. The approximations are expected to be poor, how- 
ever, in the near-tip region. 

For a plate containing a central through-thickness crack of length 2a, the 
strain at any point and the relative displacement between any two points in the 
plate may be determined by using the Westergaard method and the appropriate 
stress functions [24], The Westergaard method and the derivation of the rele- 
vant strain and displacement equations are summarized in the Appendix. The 
longitudinal strain, ^y, and displacement, v, produced by ^ remote loading on a 

finite width plate are given by Equations 9 and 10: 
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For the present discussions, these equations are assumed to be valid also for 
compression loading, since the crack is assumed to be notch-like and to remain 
separated. The longitudinal strain and displacement produced by the distributed 
closure forces are obtained by integration of the solutions for two pairs of unit 
splitting forces applied to the crack faces, Figure 14. The longitudinal strain 
and displacement by the unit splitting forces are given by Equations 11 and 12: 
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The longitudinal strain and displacement produced by the distributed crack 
closure forces are then given by Equation 13 and 14: 
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The longitudinal strain and displacement produced by the combination of 
remote loading and crack closure forces are then obtained by the superposition 
of the solutions given by Equation 15 and 16: 

Ev = ( Ev ) remo t e + ( Ev )closure ( 15 ) 

E ^y“ ( E 9 remote* ^ E 9 closure 

The following definitions of terms were used for the previous equations: 
z = x + i y ; x,y, coordinate of the point 
E, V = Young's modulus and Poisson ratio respectively 
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B, W 
b 

a close 
f(b) 

Vsec ( na/W) 


thickness and width of the specimen, respectively 

the distance between the center of the crack and splitting force 

the distance between the center of the crack and the contact 
line 

the magnitude and distribution of closure force 
correction factor for the finite width specimen 


Unfortunately, the magnitude and distribution of the crack closure forces 
are not known. Under the assumptions of the model, an estimate may be made 
on the basis of the crack closure data and the load- strain results. Before pro- 
ceeding with this estimate, qualitative features of this model may be examined 
and discussed in terms of the load-strain results. 

B. Behavior of Load-Strain Curves 

To facilitate the discussion of the load-strain curves, the longitudinal 
strains at fixed gage locations produced by unit splitting forces were computed 
as a function of the location of these forces, and are shown in Figures 15 and 16. 
Two different crack lengths were chosen. For one, all of the gages were 
located ahead of the crack tip; and the other, all behind the crack tip. The 
crack tips were sufficiently far away from the gages in both cases such that 
the approximations would be reasonable. Some specific cases are discussed 
below. 

1. Strain Gage Located Ahead of the Crack Tip 

A typical load-strain record for a strain gage located ahead of the 
crack tip is shown in Figure 17. This curve deviated from the load-strain 
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curve for a notch-like elastic crack whenever the crack closure occurred. 
Whenever there is crack closure, compressive forces acting on the crack sur- 
faces are expected. These closure forces produce tensile strain at the point 
ahead of the crack tip, just like a wedge force acting on the crack (Figure 15). 
The tensile strain produced by the closure forces makes the load-strain curves 
deviate from the idealized one. Upon further unloading, both the closed area 
and the closure force increase, so that the strain-load curve deviates further 
away from the idealized curve. The shadowed area depicted in Figure 17 repre- 
sents the effect of the closure force. 

2. Strain Gage Located Just Behind the Crack Tip 

A typical load-strain curve for a strain gage located just behind the 
crack tip is shown in Figure 18. Initial deviation from the curve for the notch- 
like elastic crack before the onset of crack closure may be attributed to the 
onset of reversed yielding [30], Once closure occurred, the curve deviated 
substantially from the idealized one. Now, the closure forces work like dis- 
tributed forces acting on a straight boundary, and produce compressive strain 
on the gage points. The shadowed area in Figure 18 depicts the influence of the 
closure force and the cross hatched area shows approximately the influence of 
reversed yielding. Since the strain caused by closure forces and by reversed 
yielding are compressive instead of tensile, the shadowed area in Figure 18 is 
on the other side of the idealized curve when compared with the shadowed area 
in Figure 17. 

It should be noted that the rate of decrease of strain was very high just 
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after the onset of crack closure, and slowed down and approached a constant with 
further unloading. The initial rapid decrease was caused by both the rapid 
increase in closed area and of the closure force effect. As the contact line 
moved further away from the strain gage, the closure force on the newly closed 
area produced little additional strain on the strain gage. Hence, the rate of 
decrease became principally a function of the remotely applied load. 

3. Strain Gage Located Far Behind the Crack Tip 

For constant amplitude loading, the electrical potential measurements 
indicated that the onset of crack closure occurred at increasingly higher loads 
with increasing K max , that is, with increasing crack lengths, Figure 12b. 

Strain gage measurements (Figure lib) appear to suggest, however, that the 
onset of crack closure occurred at progressively lower loads for the longer 
crack lengths. This apparent discrepancy was caused by the reduction in sensi- 
tivity when the strain gage is far away from the crack tip and can be explained 
in terms of the influence line depicted in Figure 16. Although the onset of crack 
closure actually occurred at a higher load for a longer crack length (as indicated 
by the electrical potential method), the strain produced at the gage location 
(when the crack tip was well past the gage) was small because of its remoteness 
from the closure forces near the crack tip. (Of course, the strain produced at 
the gage location by reversed yielding at the crack tip would have been small 
also). The strain gage could sense the effect of crack closure only when the 
contact line approached the gage. This occurred at loads much lower than that 
at the onset of closure, as shown in Figure lib. Examination of Figure 11 shows 

that the strain produced by the splitting force changed signs depending on the 
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location of the force in relation to the strain gage; being tensile over certain 
regions far removed from the strain gage and compressive near the gage. This 
change in sign, in combination with changes in magnitude with location, is con- 
sistant with the observed behavior. 

C. Estimate of Closure Forces 

An estimate of the magnitude and distribution of the closure force may be 
made on the basis of the superposition model, and experimental data on closure 
and on load-strain variations. The estimate was made for the case of R = 0. 05, 
since electrical potential measurements of crack closure and companion load- 
strain data were available. The extent of crack closure (through-thickness 
average) was determined from the electrical potential records, Figure 7(b). 
Closure forces with assumed magnitude and assumed distribution over the 
closed region were then used, in conjunction with the corresponding remote 
loading, to calculate the longitudinal strain (by numerical integration of Equation 
16) at the gage location and compared with the experimental load-strain data. 
Typically, uniform closure forces (producing yield level stresses) were applied 
to a region near the crack tip, and rapidly decaying forces were applied to the 
remainder of the closed region. This distribution was considered to be consistent 
with the physical model. The forces that produced the best fit with the experi- 
mental data (Figure 19) were considered to be a reasonable estimate of the 
closure forces. Closure forces determined by this iterative procedure are 
shown in Table 3 and appear to be reasonable. A portion of the closure force at 
the crack tip may be regarded as a contribution from the compressive yield 
zone, although such an identification was not made explicitly. 
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In arriving at these estimates, the distribution of closure forces was 
assumed to be constant in the thickness direction, and the crack front and the 
closure line were assumed to be straight and perpendicular to the specimen sur- 
face. Actually, the contact area would resemble that illustrated in Figure 20. 
This assumed three-dimensional distribution is quite reasonable, since the 
permanent tensile deformation near the surface is much larger than that in the 
center, because of the difference of plastic zone size (Figure 21). In this case, 
the distribution of the closure force must be solved as a three-dimensional pro- 
blem. Because of this two dimensional approximation, the procedure for esti- 
mating the closure forces may be used only when the strain gage is located 
close to the crack tip. 

D. Computation of Crack Opening Displacement 

The relative displacement between several pairs of points in the longitu- 
dinal direction (that is, crack opening displacement) were computed by using 
the superposition model, and the magnitude and distribution of closure force 
estimated previously. The equations used in this computation are Equations 9, 
11, 13, and 15. 

The computed displacements are depicted in Figure 22. These curves are 
similar to those observed by Elber [10], Figure 1, and suggest that the onset of 
crack closure detected by several experimental techniques would be the same, 
Figure 23. 
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V. GENERAL DISCUSSIONS 


It has been shown that the crack tip stress intensity factor, K, defined by 
linear elasticity, is the most appropriate parameter for characterizing crack 
growth [2, 3,4]. Because of a lack of understanding of the mechanism for 
fatigue crack growth, no fundamental fatigue crack growth laws have been pro- 
posed. The various "so-called" laws are empirical representations of avail- 
able experimental data, hi proposing the relationship given by Equation 1, 

Paris assumed that AK was the most significant parameter, and that stress 
ratio, R, (or, or is of only secondary importance.* In practi- 

cal cases, the effect of R can be important. Attempts have been made to incor- 
porate the influence of R into the growth rate relations, again through empiri- 
cal correlation [5, 6], 

With the observation of crack closure, it is thought that a viable physical 
basis for explaining the influence of R on fatigue crack growth had been devel- 
oped [10]. Because of closure, the effective stress intensity range, AK e £j, 
may differ considerably from the applied AK, and is a function of R. Elber 
showed that good agreement with experimental data was obtained on this basis. 
Unfortunately, several questions remain u nr esolved. Firstly, Elber argued 
that as soon as crack closure begins, the singularity at the crack tip is removed. 
Thus, the corresponding K value at the onset of crack closure is to be regarded 
as the K m i n in the load cycle; that is, a truncation of the lower portions of the 
load cycles. Experimental results from this investigation show the strains ahead of the 
* Environmental effects will not be considered in this discussion. 
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crack tip continue to decrease with unloading beyond this point, Therefore, since 
fatigue damage is related to the cyclic strain range, the effective K m j n is likely 
to be somewhat lower than the level at the onset of closure (Figure 24). These 
differences in viewpoint cannot be resolved by the approximate two-dimensional 
model developed in this report. Secondly, because of the preliminary nature of 
Elber's closure data (see Introduction), the correlation between the rate of 
fatigue crack growth and AKgff given by Equation 5 requires additional verifi- 
cation . Recent data obtained by FitzGerald and Wei [31] indicate that there is 
an effect of R on fatigue crack growth, for R values from 0. 05 to 0. 9, in the 
Ti-6A1-4V alloy used in this study, Figure 25. The present results indicate 
the absence of closure at R greater than 0.3. Hence, crack closure is likely 
to be only one of several factors that contribute to the observed R effect. 

The importance of delay, or retardation in the rate of fatigue crack growth, 
produced by load interactions in variable-amplitude loading on the accurate pre- 
diction of fatigue lives has been well recognized for some time [32,33,34]. 

/ V - 

Jonas and Wei [35] showed that the phenomenon of delay is very complex and 
can depend on a broad range of loading variables. Due to a lack of physical 
and phenomenological understanding of the effects of load interactions on fatigue 
crack growth, no successful model has been proposed to account for delay. 

In general, delay effects have been attributed to such things as crack tip 
blunting [36], residual stress ahead of the crack tip [34], and more recently, 
crack closure proposed by Elber [10]. The first two of these theories do not 
lend themselves to direct experimental observation. The third, crack closure 
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proposed by Elber, may be verified experimentally by determining if delay 
occurs in the absence of closure. For a loading spectrum depicted on Figure 
26, there should not be any crack closure occurring for K m ^ n equal to or 
greater than 12 ksi -/in, (R > 0. 3), according to the result of the last section. 
However, significant amount of delay was observed, Figure 26. This 
experimental result shows again that crack closure cannot be regarded as the 
sole cause for delay. 

The present series of experiments have shown that crack closure can and 
does occur. Its occurrence is dramatically evidenced by changes in electrical 
potential, strain and crack opening displacement with load. However, crack 
closure cannot account for all of the observed fatigue behavior. Additional 
research is needed to better establish its significance. 
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VI. CONCLUSION 


The following conclusions can be made on the basis of the experimental 
results on the Ti-6A1-4V alloy: 

1. Crackclosure does occur during fatigue. For the Ti-6A1-4V alloy, 
closure was observed at stress ratios, R, between 0.05 and 0.3. 

2. The extent of closure can be determined by electrical-potential mea- 
surement method with reasonable accuracy. 

3. The stress intensity factor at the onset of closure depends on the stress 
ratio, R, and on K max . 

4. Onset of closure can be detected by measurements of such quantities 
as electrical potential, crack-opening displacement, and strain. 

5. The closure force near the crack tip was estimated to produce yield 
stresses. A part of this estimated closure force may be attributed to the 
plastic zone. 

6. Crack closure cannot be regarded as the sole cause for the various 
observed phenomena for fatigue. For example, it cannot completely account 
for the influence of stress ratio, R, on fatigue crack growth; neither can it 
fully explain the delay phenomenon under variable- amplitude loading. 
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TABLE 1 


Chemical Composition and Tensile Properties 
(Reactive Metal Ingot No. 293831, Lot 05) 

Chemical Composition - Weight Percent 
(Ingot Analysis) 


C N 

Fe A1 V 

O H (ppm) 

Ti 

\ 

0.03 0.013 

0.13 6.2 4.2 

0. 12 90 

Balance 


Tensile Properties* 

Tensile 

Yield Strength Strength 

Elongation 
in 2 in. 

Direction 

ksi (MN/m2) 

ksi (MN/m2) 

percent 

Transverse 


149.5 (1,031) 

152. 6 (1,052) 

12. 5 


149.8 (1,033) 

153.4 (1, 058) 

12. 5 


148.3 (1,023) 

151.4 (1,044) 

11. 7 

(Average) 

149. 2 (1, 029) 

152. 5 (1,051) 

12. 2 

Longitudinal 


— 

150. 7 (1, 039) 

11. 7 


141. 1 (973) 

151. 0 (1,041) 

11. 7 


141.9 (978) 

151.8 (1,047) 

12.5 

(Average) 

141.5 (976) 

151.2 (1,042) 

12.0 


* Production annealed 1,450° F, 15 min. + air cool. 
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Figure 4. Schematic Diagram of Electrical- Potential System for Crack Measurement 
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Figure 5. Schematic Diagram of Environmental Control System 
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Figure 6a. Load versus Electrical Potential Curves (before 
the Compressive Loading Cycle) (R = 0. 1) 
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Figure 6b. Load versus Electrical Potential Curves (before 
the Compressive Loading Cycle) (R = 0. 2) 
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Figure 6c. Load versus Electrical Potential Curves (before 
the Compressive Loading Cycle) (R = 0. 3) 
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Figure 6d. Load versus Electrical Potential Curves (before 
the Compressive Loading Cycle) (R = 0. 5) 
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Figure 8. Load versus Electrical Potential Curves (after the 
Compressive Loading Cycle) 
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Figure 11a. Load-strain Curves for Gage A. 

(x = 0. 347 in; y = 0. 036 in. , see Table 2) 
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Figure lie. Load-Strain Curves for Gage D. 

(x = 0. 797 in. ; y = 0. 079 in. , see Table 2) 




58 


Figure 12a. Relationship between U and K for Different Stress Ratios 
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Figure 13. Superposition Model 
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Figure 16. Relationships between Strains and the Position of Unit 
Splitting Forces with Strain Gage Located Behind the 
Crack Tip (a = 0. 908 in. ) 
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Figure 17. Typical Load versus Strain Curve for a Strain Gage 
Located Ahead of the Crack Tip 



Figure 18. Typical Load versus Strain Curve for a Strain Gage 
Located Just Behind the Crack Tip 
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Figure 19. Longitudinal Strain Calculated by the Superposition 
Model and the Estimated Closure Force 



Closure Line 

( Shaded Area = Aa x b ) 


Figure 20. Comparison of Probable Closure Line with Closure Line 
Indicated by the Electrical Potential Measurements 
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Figure 22. Relationship between Load and Relative Displacement 
between Gage Points Computed by Using the Super- 
position Model 
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C.O.D. MEASUREMENT (Elber) STRAIN MEASUREMENT 

Figure 23. Schematic Shown of the Onset of Crack Closure Detected 
by Three Different Experimental Techniques 
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APPENDIX 


The equation of equilibrium for two dimensional elastic problems can be 
found in any elementary text of the theory of elasticity. 

- 0 


_36x _&£y _ 

dX dj) 


d d<$y _ 


- 0 


A-l 


ax ' dj 
Txy - Tyx 

The strain-displacement relationships and Hook's law lead to the compatibility 
equation. 

3* 


(a? + d?) ■ 0 


A- 2 


As is easily checked, equations A-l are satisfied by an Airy's stress function, 
$ (x, y), if the stress components are defined as follows: 

<5x = 

^ = dX z 


d z & 

Zf 

d 1 ! 


A-3 


7xy — ~ 




Substitutions of Equation A-3 into Equation A- 2 gives the biharmonic equation 
in 3E . 


V 4 ± = 


dtf 


0 , #I_ 

3X z dP dp 


= o 


A-4 


For a specific problem, the stress function, £& » must satisfy Equation A-4 and 
the appropriate boundary conditons. Theory of functions of complex variables 
provides a powerful tool for the solution of two-dimensional elasticity problems. 
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Define a complex variable, z, by 

z = x + i y A-5 

If a function of this complex variable, Z(z) > is analytic then its derivatives 
exist, that is, 


z = 

dz 

7 — 

dZ 

4J — 

dz 

7' = 

dz 

Li — 

dz 


Let Z = ReZ + ilm Z, then the Cauchy Rieman conditions state that: 


A- 6 


It follows that 


dX 


ReZ = 


Jt 


Im Z = Re Z 



ImZ 




Re Z = Im Z 


V 2 (ReZ) (ImZ) = 0 


A- 7 


Similarly, all derivatives of Z have harmonic real and imaginary parts. If 
the '/l are each harmonic, that is 

V 2 )/'- = 0 A-8 

then any stress function in the form 

3 = Vl + x % + y% + ( x2 + 5^ A " 9 


will automatically satisfy the Equation A-4. 

In conformity with Equations A-5 - A-9, Westgaard defined an Airy 
stress function, I , by 
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A-10 


3i = Re Z + y Ira Z 

This function satisfies the equilibrium Equations, A-l, and the compatibility 
condition, Equation A-2. 

Using Equation A-3, A-7, and A-10 the stress are defined by Equation A-ll: 


(5 X = Re Z -y ImZ' 
= Re Z + y Im Z* 

T xy = -y ReZ ’ 


A-ll 


The displacement- strain relations are given by Equation A-12: 



3 v 




A-12 


The strain-stress relations for generalized plane stress condition are given by 
Equation A- 13 


E ~ dy ~ Vdx 


A- 13 


Combining Equation A-ll, A-12, and A-13, the in-plane displacements and 
strains are obtained. 


Eu = ( 1 - y) Re Z - ( 1 + V) y Im Z 
Ev = 2ImZ - (l+V)yReZ 


A- 14 
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A-15 


Ee x = (l-y)ReZ-(l+V)yImZ 


E€ y = ( 1 -if) Re Z + ( 1 +y> y Im Z 


For a central crack of length 2a in an infinite plate with uniform biaxial stress € 
at infinity, the stress function is given by Equations A-16 [24]: 


Z = 0 Jz 2 -a 2 


Z = 6 z//z 2 -a 2 

Z *= 6/ /z 2 - a 2 - 6z 2 / /(z 2 - a 2 )" 


A-16 


The longitudinal strain and displacement at a point A(x, y) are then obtained 
from Equations A-14 and A-15. 


Ee y = (5- [(!->>) Re(- /: ^--) + (1+V)ylm 


(z 2 - a 2 )' 


Ev = <5 [2 Imj/z 2 - a 2)_ (i +*) y Re j *■— ) 


A-17 


The stress function for a central crack with two equal pairs of splitting forces 
at x = +b (Figure 14) is given by Equations A-18. 


- _ 2 P -1 z 2 -a 2 

z ~ Tr tan fir? 


J 2 v 2 

z ia -b 


(z2-b 2 ) 


/ 2 2 
z -a 


A-18 


2 P ^J ^-b 2 , 2z 2 /a 2 -? . 

( Z 2 -^)P^T? (z 2 -^) 2 ^ 7 ^" 
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The longitudinal strain and displacement at a point A(x, y) is then given by 
Equations A- 19. 


* V 4r 

y 7t 




a^- 




(z 2 - b 2 ) j (z 2 - a 2 )^ 
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( z^- h^)^Jz^~ 
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